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The rate of period change in pulsating DB white dwarf stars 
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Abstract. In this work, we present the theoretically expected rates of pulsation period change for V777 Her (DBV) 
variable stars. To this end we employ new evolutionary models representative of pulsating DB white dwarf stars 
computed in a self-consistent way with the predictions of time-dependent element diffusion. At the hot edge of the 
DB instability strip, the envelopes of the models are characterized by a diffusion-induced double-layered chemical 
structure. We compute the numerical values of rates of period change by solving the equations of linear, adiabatic, 
nonradial stellar oscillations. We examine the effects of varying the stellar mass, the mass of the helium envelope 
and the neutrino emission on the expected period changes. We present extensive tabulations of our results which 
could be useful for comparison with future detections of the rate of period change in pulsating DB white dwarfs. 
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1. Introduction 

V777 Her (or DBV) stars are (;(gravity)-mode variable 
white dwarfs with helium-rich atmospheres (DB) and in- 
termediate effective temperature (Tcff ^ 25000 K), the 
pulsating nature of which was predicted two decades ago 
on theoretical grounds by Winget et al. (1982a) and 
shortly after confirmed observationally by Winget et al. 
(1982b). Since then, considerable effort has been devoted 
to studying these stars. In particular, the multiperiodic 
star GD 358, the most extensively studied member of the 
DBV class, has been the subject of numerous investiga- 
tions devoted to disentangling its internal structure and 
evolution, initially by means of "hand on" asteroseismo- 
logical procedures (Bradley & Winget 1994) and later by 
employing objective fitting techniques (see, e.g, Metcalfe 
et al. 2000; 2001; 2002). In particular, Metcalfe et al. 
(2001) — see also Metcalfe et al. (2002) — have recently 
applied genetic algorithm-based procedures to place con- 
straints on the ^^C(q;,7)^^0 reaction rate from inferences 
for the abundance of central oxygen in GD 358. 



As a variable white dwarf cools down, its oscillation 
periods (P) vary in response to evolutionary changes in 
the mechanical structure of the star. Specifically, as the 
temperature in the core of a white dwarf decreases, the 
plasma increases its degree of degeneracy and the Brunt- 
Vaisala frequency — the most important physical quantity 
in g-mode pulsations — diminishes, and the pulsational 
spectrum of the star shifts to longer periods. On the other 
hand, residual gravitational contraction (if present) acts 
in the opposite direction, thus shortening the pulsation 
periods. Competition between the increasingly degeneracy 
and gravitational contraction gives rise to a detectable 
temporal rate of change of periods (P = dP/ dt) . Roughly, 
the rate of change of the pulsation period is related to the 
rates of change of the temperature at the region of the 
period formation, T, and of the stellar radius, i?*, by the 
expression (Winget et al. 1983): 



P 
P 



-a h — , 
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(1) 
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where a and b are constants whose values depend on the 
details of the white dwarf modeling {a,b « 1). The first 
term in Eq. corresponds to the rate of change in pe- 
riod induced by the cooling of the white dwarf and it 
is a positive contribution, while the second term repre- 
sents the rate of change due to gravitational contraction 
and it is a negative contribution. In principle, the rate of 
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TABLE 1 



Periods and Relative Rates of Period Change (d — 1 modes) for selected models with M* — 0.50Mq and Mhc = 8 x 10 ^Af« 
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1.88 


585.02 


2.02 


599.78 


1.91 


12 
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13 
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14 
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15 
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4.88 
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796.58 


1.84 


811.81 


1.84 


830.30 


1.72 


17 


782 


04 


6.04 


798.52 


4.06 


808.98 


1.88 


817.49 


1.94 


834.90 


2.71 


858.10 


2.51 


881.87 


1.94 


18 
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20 
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21 
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NOTE - All P/P values are in units of 10"^^ s"\ 



change of period can be measured by observing a pulsating 
white dwarf over several seasons when one or more stable 
pulsation periods are present in their light curves. In the 
case of pulsating DA (hydrogen-rich atmospheres) white 
dwarfs — also termed DAV or ZZ Ceti variable stars — 
cooling dominates over gravitational contraction, in such 
a way that the second term in Eq. is negligible, and 
only positive values of the rate of change of period are ex- 
pected. This is the case of the ZZ Ceti G117-B15A, for 
which Kepler et al. (2000) have quoted a rate of change 
of the 215.2 s period of P = (2.3 ± 1.4) x lO'^^ s s'^ In 
addition, Mukadam et al. (2003) report that the period at 
213.13 s observed in other ZZ Ceti star, R548, drifts at a 
rate P < (5.5 ± 1.9) x 10"^^ s s~^ For the high effective 
temperatures characterizing the DOV/PNNV instability 
strip, gravitational contraction is still significant, to such 
a degree that its influence on P can overcome the effects 
of the cooling. In this case the second term in Eq. is 
not negligible anymore and consequently positive or nega- 
tive P values are possible. For these hot pre-white dwarfs, 
what determines that P is either positive or negative is 
the character of the mode: if their pulsation period corre- 
sponds to a confined mode^ then a positive value of P is 
expected, while if the mode is a trapped one^ we would 
expect a negative P value (see Kawaler & Bradley 1994). 
For PG1159-035, the prototype star of the DOV class, a 
rate of P = (13.07±0.03) x 10"" s s'^ has been measured 

^ That is, a mode with their eigenfunctions sensitive to the 
deep regions of the star. 

^ That is, a mode whose eigenfunctions are concentrated 
mostly toward the outer layers. 



for the 516 s pulsation period, although uncertainties in 
the modeling of its interior structure have led to ambigu- 
ous conclusions (Winget et al. 1985; Costa et al. 1999). As 
for the DBV stars, in which the influence of the gravita- 
tional contraction on P is negligible, the expected rate of 
period change from theoretical evolutionary computations 
should be positive and of the order of ~ 10^^^ — 10^^^ s 
s"""^, although detections of P have not yet been made in 
any DBV star. 

Observational measurement of P provides a sensitive 
probe of the structure and evolution of white dwarf stars. 
As shown by Winget et al. (1983), a measurement of the 
rate of period change of a pulsating white dwarf consti- 
tutes, particularly for DAV and DBV white dwarfs (which 
evolve at almost constant radius) a direct measurement 
of the cooling rate of the star. This, in turn, provides 
valuable information about the core chemical composi- 
tion. Also, measurement of the rate of period change of 
pulsating white dwarfs in the DAV, DBV and DOV in- 
stability strips would allows astronomers to calibrate the 
cooling sequence age. This, in turn, could be employed to 
infer the age of the galactic disk in the solar neighbor- 
hood (Winget et al. 1987). In addition, the measurement 
of P in variable white dwarfs can be employed to set con- 
straints on particle physics. McGraw et al. (1979) were the 
first to suggest that hot pulsating white dwarfs could be 
employed to determine the effect of neutrino cooling as a 
star becomes a white dwarf (see also Winget et al. 1983). 
The influence of neutrino energy loss on P was discussed 
in detail by Kawaler et al. (1986) for the case of DBV 
and DOV white dwarfs. In addition, Isern et al. (1992) 
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TABLE 2 



Periods and Relative Rates of Period Change {£ — 1 modes) for selected models with Af, — 0.60Mq and A/hc = 8 x 10 ^Af* 
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4.10 
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1.00 
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569.46 


4.28 


581.40 


3.30 


592.27 


2.41 
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1.71 


615.98 


1.41 
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1.20 
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3.52 
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2.40 


627.25 


2.06 


640.25 


1.99 


655.81 


1.75 
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NOTE - AU PjP values are in units of 10"^'' 8"^ 



— see also Corsico et al. (2001a) — have explored the ef- 
fect of axion emissivity in DAV stars. Recently, O'Brien & 
Kawaler (2000) have discussed the possibility of inferring 
limits on the theoretically determined plasmon neutrino 
emission rates by employing DOV white dwarfs. 

More recently, D. Winget and collaborators have 

drawn attention the attractive possibility of employing 
DBV stars as plasmon neutrino detectors (Winget et al. 
2004). These authors — see also Kawaler et al. (1986) — 
have noted that in the hotter region of the DBV instability 
strip, the plasmon neutrino energy losses are several times 
larger than the losses due to photon emission. Their re- 
sults suggest that measurement of the evolutionary period 
change in hot DB white dwarf stars would constitute an 
excellent probe of the plasmon neutrino production rates, 
if period changes in DBV stars could be assessed. The au- 
thors discuss several observational strategies to estimate 
P for DBV white dwarfs, in addition to the ongoing ob- 
servations of the hot DBV EC 20058 by Sullivan et al. 
(2004). 

It is worth noting that the study of Winget et al. (2004) 
is based on stellar models for pulsating DB stars with a 
chemical structure characterized by a pure helium enve- 
lope atop a carbon-oxygen core, i.e., a single- layered en- 
velope structure. However, on the basis of evolutionary 
calculations including time-dependent clement diffusion, 
Dehner & Kawaler (1995) and Gautschy & Althaus (2002) 
have found that, if DB white dwarfs descend from PG1159 



stars^, their envelopes would be characterized by the pres- 
ence of a double-layered chemical structure. Indeed, such 
calculations show that by the time the DB instability 
strip is reached, models are characterized by two differ- 
ent chemical transition zones, i.e., a double-layered con- 
figuration. In fact, above the carbon-oxygen core, there 
exists an envelope consisting of an intershell region rich 
in helium, c;arboii and oxygen, the relics of the short-lived 
mixing episode occurred during the last helium thermal 
pulse that leads to the born-again episode, and an over- 
lying pure helium mantle resulting from the gravitational 
settling of carbon and oxygen. More recently, Fontaine 
& Brassard (2002) have demonstrated that the theoreti- 
cal period spectnun of DBV white dwarf models which 
incorporate a double-layered envelope, turns out to be 
markedly distinct from that expected for a single-layered 
configuration. As shown by these authors, this is partic- 
ularly important when attempts at constraining the core 
composition of DBs and the ^^C(a, 7)^^0 reaction rate are 
made. This point has recently been addressed by Metcalfe 
et al. (2003), who have incorporated both the double- 
layered envelope feature and adjustable carbon-oxygen 
cores in DB asteroseismological fittings. 



^ PG1159 stars arc hydrogen-deficient post-asymptotic giant 
branch (post-AGB) stars that are widely believed to be the 
result of a born-again episode (see Schonberner 1979 and Iben 
et al. 1983) and the immediate predecessors of most hydrogen- 
deficient white dwarfs. 
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TABLE 3 



Periods and Relative Rates of Period Change {£ — 1 modes) for selected models with Af, — 0.60Mq and A/hc = 9 x 10 '^Mt 
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311.14 
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325.56 


3.28 
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1.82 
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349.21 


3.97 


355.58 


3.25 


361.91 
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1.61 
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1.26 
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576.74 


4.63 


587.92 


3.21 
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629.52 
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1.16 
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4.09 


615.16 


3.41 


627.28 


3.02 


643.50 


2.60 


660.97 


2.09 


678.25 


1.52 


698.82 


1.02 
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638.71 


3.32 


648.61 


2.88 


659.84 


2.72 


675.37 


2.38 


692.07 


1.89 


708.57 


1.42 


731.31 


1.19 


16 


675.30 


3.89 


688.87 


3.78 


703.85 


3.24 


721.53 


2.32 


736.43 


1.42 


748.60 


0.99 


768.41 


1.12 


17 


715.66 


5.05 


732.33 


4.10 


748.03 


3.00 


763.80 


1.86 


777.88 


1.47 


794.60 


1.40 


821.62 


1.28 


18 


753.44 


4.67 


768.98 


3.53 


782.79 


2.52 


798.90 


2.12 


818.05 


1.95 


839.09 


1.55 


867.38 


1.25 


19 


786.84 


4.59 


801.97 


3.21 


815.80 


2.61 


834.84 


2.43 


856.35 


1.98 


877.28 


1.41 


905.36 


1.27 


20 


821.97 


2.86 


833.66 


2.84 


849.31 


3.05 


871.74 


2.62 


894.04 


1.86 


913.54 


1.26 


943.14 


1.34 


21 


861.94 


3.76 


878.81 


3.69 


897.11 


3.06 


917.38 


2.02 


933.75 


1.29 


951.07 


1.27 


984.13 


1.40 


22 


900.86 


4.99 


922.07 


4.19 


942.30 


3.05 


962.10 


1.85 


980.79 


1.61 


1002.96 


1.41 


1036.37 


1.35 


23 


932.29 


4.99 


953.51 


3.96 


972.99 


2.88 


995.86 


2.42 


1022.65 


2.13 


1050.03 


1.55 


1085.76 


1.38 


24 


975.02 


4.44 


992.10 


2.77 


1006.48 


2.25 


1028.72 


2.40 


1055.00 


1.95 


1080.16 


1.39 


1118.25 


1.44 



NOTE - AU PjP values are in units of 10" fT^ . 



In a recent work, Althaus & Corsico (2004) have pre- 
sented new evolutionary models of DB white dwarf stars 
for various masses of the helium content and several stel- 
lar masses, computed in a self-consistent way with the 
predictions of time-dependent element diffusion. The ini- 
tial outer layer chemical stratification assumed for such 
models corresponds to that characterizing PG1159 stars. 
By the time the domain of the DBVs is reached, the en- 
velopes of these models are characterized by a double- 
layered chemical structure induced by diffusion. The au- 
thors find that, depending on the stellar mass, if DB 
white dwarf progenitors are formed with a helium con- 
tent smaller than ~ 10~^M*, a single-layered configura- 
tion is expected to emerge during the DB pulsation in- 
stability strip. As shown in that paper, the period spac- 
ing diagrams exhibit mode-trapping substructures when 
a double-layered configuration characterizes the envelope 
of the models, substructures that are virtually absent in 
single-layered envelope models. 

In view of the Winget et al. (2004) 's claims about the 
potential of employing DBV white dwarf stars to place 
constraints on the plasmon neutrino emissivity, we judge 
that the computation of the rate of period change in the 
frame of new stellar models which incorporate an updated 
input physics, particularly time-dependent diffusion pro- 
cesses as well as realistic initial envelope chemical strat- 
ifications as predicted by the evolutionary history of the 
progenitor stars, would be worthwhile to be done. This 
is the aim of the present paper. An additional motiva- 



tion is the lack of modern tabulations of the rate of pe- 
riod change for DBV stars in the literature. Specifically, 
we present theoretical values for the rates of pulsation 
period change for the effective temperature range of in- 
terest. This is done by means of several tables providing 
P and PjP (= dlnP/dt) values corresponding to dipole 
modes {£ = 1) for different stellar masses, effective tem- 
peratures and thickness of the helium envelope. DB white 
dwarf evolutionary models are the same as those presented 
in Althaus & Corsico (2004). We compute the numerical 
values of P/P by solving the equations of linear, adia- 
batic, nonradial stellar oscillation for evolutionary models 
representing DB white dwarfs. We briefly examine the ef- 
fect of changing the various structural parameters on the 
theoretical rate of period change. 

The rest of the paper is organized as follow. In §2 we 
describe our evolutionary DB white dwarf models, and in 
§ 3 we examine the effects of varying the stellar mass, the 
mass of the helium envelope and the neutrino emission on 
the P/P values. Finally, we close the paper with a short 
summary in §4. 

2. Input physics and evolutionary models 

The evolutionary models employed in this work have been 
obtained with the DB white dwarf evolutionary code de- 
veloped at La Plata Observatory. The code is that de- 
scribed in Althaus & Corsico (2004) and references therein 
(see also Gautschy & Althaus 2002). In particular, mi- 



Corsico & Althaus: The rate of period change in DBV white dwarf stars 5 

TABLE 4 



Periods and Relative Rates of Period Change {£ — 1 modes) for selected models with Af, — 0.60Mq and Muc = 1 x 10 ''M* 









27956 K 


26901 K 


25907 K 


24969 K 


24086 K 


23095 K 


22021 K 


k 


P 


s] 


P/P 


P \a] 


P/P 


P fsl 


PIP 


P fsl 


PIP 


P fsl 


PIP 


P H 


PIP 


P fsl 


PIP 


1 


128 


01 


6.60 


131.90 


5.43 


135.85 


4.42 


139.85 


3.56 


143.90 


2.84 


148.82 


2.18 


154.65 


1.62 


2 


165 


08 


4.84 


168.74 


3.99 


172.45 


3.27 


176.21 


2.65 


180.02 


2.14 


184.64 


1.64 


190.08 


1.23 


3 


196 


61 


4.98 


201.14 


4.18 


205.84 


3.48 


210.68 


2.87 


215.66 


2.34 


221.74 


1.80 


228.90 


1.33 


4 


235 


43 


5.20 


241.07 


4.32 


246.82 


3.54 


252.65 


2.86 


258.49 


2.27 


265.43 


1.69 


273.34 


1.21 


5 


276 


42 


5.28 


283.08 


4.31 


289.75 


3.47 


296.38 


2.75 


302.94 


2.17 


310.74 


1.63 


319.88 


1.22 


6 


310 


29 


4.75 


316.96 


3.86 


323.68 


3.14 


330.44 


2.53 


337.28 


2.05 


345.66 


1.60 


355.69 


1.21 


7 


350 


48 


4.94 


358.32 


3.98 


366.00 


3.14 


373.48 


2.44 


380.66 


1.86 


388.78 


1.32 


397.58 


0.92 


8 


392 


79 


4.81 


401.05 


3.68 


408.79 


2.77 


416.04 


2.10 


423.03 


1.65 


431.50 


1.31 


442.12 


1.06 


9 


424 


64 


3.58 


431.67 


3.04 


439.23 


2.70 


447.61 


2.41 


456.84 


2.10 


468.89 


1.73 


483.91 


1.35 


10 


456 


90 


4.14 


465.85 


3.61 


475.36 


3.08 


485.36 


2.58 


495.77 


2.14 


508.55 


1.65 


523.18 


1.17 


11 


492 


46 


4.54 


503.00 


3.93 


514.11 


3.31 


525.56 


2.71 


536.96 


2.11 


549.83 


1.46 


563.35 


1.01 


12 


532 


52 


4.74 


544.15 


3.93 


555.70 


3.08 


566.51 


2.26 


576.14 


1.60 


586.75 


1.19 


600.74 


1.10 


13 


575 


36 


5.07 


587.85 


3.70 


598.60 


2.51 


607.87 


1.83 


617.27 


1.61 


630.47 


1.46 


648.41 


1.26 


14 


605 


50 


4.33 


616.96 


3.35 


628.59 


2.90 


641.70 


2.65 


656.35 


2.31 


674.87 


1.80 


696.15 


1.32 


15 


640 


04 


3.14 


649.90 


2.97 


661.76 


2.89 


675.45 


2.59 


689.99 


2.13 


706.97 


1.51 


726.02 


1.22 


16 


678 


49 


4.30 


692.81 


3.95 


708.19 


3.31 


723.29 


2.49 


736.56 


1.69 


749.71 


1.10 


767.53 


1.24 


17 


716 


50 


5.11 


733.51 


4.26 


750.04 


3.20 


764.58 


2.19 


777.50 


1.67 


794.51 


1.50 


818.67 


1.45 


18 


755 


05 


5.10 


771.67 


3.75 


785.91 


2.55 


799.06 


2.11 


814.39 


2.04 


835.69 


1.71 


862.65 


1.48 


19 


788 


57 


4.54 


803.40 


3.20 


817.64 


2.75 


834.31 


2.62 


853.05 


2.25 


875.60 


1.63 


901.48 


1.42 


20 


825 


12 


3.06 


837.40 


2.94 


853.36 


3.07 


872.03 


2.72 


890.97 


2.08 


911.30 


1.37 


937.25 


1.49 


21 


865 


39 


3.91 


882.86 


3.86 


902.12 


3.23 


920.33 


2.29 


935.32 


1.49 


952.37 


1.28 


980.51 


1.56 


22 


901 


73 


5.28 


924.05 


4.45 


945.61 


3.27 


963.81 


2.14 


980.44 


1.78 


1003.48 


1.57 


1035.34 


1.58 


23 


934 


95 


5.08 


956.02 


3.91 


974.90 


2.78 


993.43 


2.45 


1015.47 


2.31 


1044.50 


1.82 


1079.78 


1.60 


24 


977 


94 


4.74 


996.39 


3.03 


1012.01 


2.40 


1030.75 


2.45 


1052.67 


2.14 


1078.27 


1.47 


1110.15 


1.59 



NOTE - AU PjP values are in units of 10"^'' s'^ 



crophysics includes an updated version of the equation 
of state of Magni & Mazzitelli (1979), OPAL radiative 
opacities for arbitrary metallicity (Iglesias & Rogers 1996) 
including carbon- and oxygen-rich compositions, and up- 
to-date neutrino emission rates and conductive opacities. 
In particular, opacities for various metallicities are re- 
quired because of the metallicity gradient that develops 
in the envelopes as a result of gravitational settling. In 
this work, convection is treated in the framework of the 
mixing length theory as given by the ML2 parameteriza- 
tion (sec Tassoul et al. 1990). Our evolutionary models are 
calculated self-consistently with the predictions of time- 
dependent element diffusion. Our treatment of diffusion 
for multi-component gases includes gravitational settling 
and chemical and thermal diffusion for the nuclear species 
^He, ^^C and ^^O. As for the chemical composition of the 
core, we have adopted the chemical profiles of Safaris et 
al. (1997). 

Our starting stellar configurations correspond to hot 
white dwarf structures with a realistic outer layer chemi- 
cal stratification appropriate to that of hydrogen- deficient 
PG 1159 stars. These stars, presumed to be the direct an- 
cestors of most DB white dwarfs, are widely believed to be 
the result of a born-again episode experienced by a post- 
AGB remnant on the early white dwarf cooling branch. As 
a result of such an episode, most of the hydrogen content 
is completely burnt, and the envelope is eventually charac- 
terized by an uniform chemical composition of helium, car- 
bon and oxygen. In our models we have assumed, for such 



compositions, mass fractions of 0.42, 0.36 and 0.22, respec- 
tively, following observed abundance patterns in PG1159 
stars. Also, we have varied the stellar mass in the range 
0.50 — 0.85 Mq, and the mass of the helium content to be 
Mne: 8 x lO'^, 9 x lO"* and 1 x lO"* M*. With regard to 
the helium envelope, some words are in order. The mass of 
the helium envelope is constrained by the theory of post- 
AGB evolution to be ^ IQ-^M, (Iben 1989; D'Antona & 
Mazitelli 1991). Also, recent full evolutionary calculations 
predict the helium content to range from 10~^Mh< (for a 
0.6 Mqj white dwarf remnant, see Hcrwig ct al. 1999) to 
10"3M* (for a 0.93 M© white dwarf, see Althaus et al. 
2003). Such values are upper limits in the sense that post- 
AGB mass loss episodes coiild rcdiicc them considerably. 
In fact, as recently emphasized by Werner (2001), the ex- 
istence of mass-loss rates in the range 10"'' — lO~®M0/yr 
cannot be discounted in many luminous PG 1159 stars. In 
addition, tentative evidence for the persistence of mass- 
loss rates of the order 10"'' — lO"^*^ Mq/jv down to the 
domain of hot helium-rich white dwarfs has been presented 
(see Werner 2001). Thus, the existence of hot hydrogen- 
deficient pre-white dwarf stars with a helium content of 
the order of 10"'^ M* could not be discounted even in the 
case of DB white dwarf remnants with masses as low as 
0.6 MQ(see Althaus & Corsico 2004). 
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Log(l-M/MJ 



Fig. 1. The run of the squared Brunt- Vaisala frequency 
in terms of the outer mass fraction, corresponding to a 
O.6O-M0DB white dwarf model at T^s = 25338 K with 
a helium content of Mhb « 9 x 10~*M,. Dark regions 
denote the contributions of the Ledoux term B (shown in 
the inset) to the Brunt- Vaisala frequency. The chemical 
profile is displayed in the upper zone of the plot. 



3. Rates of period change 

In what follows, we describe the main results of our pul- 
sation calculations. We mostly concentrate on the theo- 
retical rate of period change. We assess the P/P values 
with the help of the pulsational code employed in Corsico 
et al. (2001b), appropriately modified to study pulsating 
DB white dwarfs. In particular, the treatment we follow to 
assess the Brunt- Vaisala frequency (N) is that proposed 
by Brassard et al. (1991). 

We begin by examining Fig. ^ in which a representa- 
tive spatial run of the Brunt- Vaisala frequency of a DB 
white dwarf is displayed. The model is characterized by 
a stellar mass of 0.60 Mq, an effective temperature of 
w 25300 K and a hchum content of Mhc « 9 x 10"'' M*. 
In addition, the plot shows the internal chemical stratifi- 
cation of the white dwarf model, and for illustrative pur- 
poses the profile of the Ledoux term B (inset), an impor- 
tant quantity related to the computation of the Brunt- 
Vaisala frequency. The figure emphasizes the role of the 
chemical interfaces on the shape of the Brunt- Vaisala fre- 
quency. In fact, each chemical transition region produces 
clear and distinctive features in N, which eventually are 
responsible for the mode trapping properties of the model. 
At the core region, the dominant feature at logg « —0.4 
((7=1 — Mr/M^,) is the result of the steep variation of the 
oxygen/carbon profile. This feature causes strong trap- 
ping of certain modes in the core region — 'core-trapped' 



modes; see Althaus et al. (2003) in the context of massive 
DA white dwarf models. Such modes are characterized by 
an unusually large oscillation kinetic energy. The less pro- 
nounced bump in N at log g « — 1 is much less relevant to 
the structure of the period spectrum. Finally, in the enve- 
lope of the model we find a double-layered chemical struc- 
ture. Despite the fact that this structure is modeled by 
diffusion — and as such, characterized by a very smooth 
functional shape — its influence on the mode trapping 
properties is by no means negligible. Indeed, Althaus & 
Corsico (2004) have found that the double-layered con- 
figuration is responsible for substructures in the period- 
spacing diagrams. Finally, the steep drop in the Brunt- 
Vaisala frequency at logg « —5.5 is caused by the opacity 
change due to the metallicity gradient induced by diffu- 
sion in the outer layers. This feature occurs close enough 
to the stellar surface and has not appreciable effects on 
the model period spectrum. 

In this work the rate of period change is estimated as 
simple differencing of the periods of successive models in 
each evolutionary sequence. We present P and P/P val- 
ues corresponding to a modest parametric survey of the 
DB evolutionary models presented in Althaus & Corsico 
(2004). In Tables 1 to 9 we provide a set of P and P/P 
theoretical values corresponding to dipole modes (^ = 1) 
for models with several stellar masses, effective tempera- 
tures and thickness of the helium envelope. We feel that 
these tables could be useful for comparison with future 
observations of the rate of period change in pulsating DB 
white dwarfs. Here, we show a sample of our rate of pe- 
riod change values. A more comprehensive tabulation is 
available at our website 

Our ^ = 1 values of P/P range from ^ 4 to ^ 7 x 10"^^ 
s-^ at Teff « 30000 K and from ~ 1.5 to - 3 x lO'^^ s'^ 
at Tcff K. 25000 K, for periods in the interval « 100-1000 s 
in the case of a representative O.6O-Af0 model, irrespective 
of the helium content. Our P/P values are almost twice 
as large as the values quoted by Kawaler et al. (1986) for 
a DB model with 0.60 Mq and Mne « 3.3 x lO'^M* at 
Teff f« 30000 K. The £ = 1,2 values of P for our mod- 
els are about 4 times greater than those of Bradley & 
Winget (1991) (their Table 12). This is mainly due to 
the distinct input physics employed and in particular to 
the fact that these authors do not consider neutrino emis- 
sion in their Tassoul et al. (1990) 's carbon-core DB white 
dwarf models. As compared with the values of Bradley 
et al. (1993) (their Table 6, corresponding to a O.6-M0 
carbon-core white dwarf model), our 1 = 2 values of P 
are about 1.55 times greater, depending on the effective 
temperature. Unfortunately, the lack of extensive tabula- 
tions of P in the Bradley et al. (1993)'s paper prevent us 
from making a comprehensive comparison between their 
results and ours. We note that our models make use of 
a much more updated description of the opacities, rates 
of neutrino emission and core chemical composition than 
that used in Bradley et al. (1993). 

^ http : //www . f caglp . unlp . edu . ar/ evolgroup/ 
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TABLE 5 

Periods and Relative Rates of Period Change {£ — 1 modes) for selected models with M* — 0.65Mq and Mhc = 9 x 10~*M« 





Tcff 




28041 K 


27090 K 


26020 K 


25024 K 


23943 K 


22945 K 


22039 K 


k 


P 


si 

J 


P/P 


P [si 


PIP 


P [si 


PIP 


P [si 


PIP 


P [si 


PjP 


P [si 


PIP 


P [si 


PIP 


1 


126 


88 


5.11 


130.38 


4.06 


134.59 


3.05 


138.82 


2.29 


143.81 


1 


69 


148.90 


1.31 


153.99 


1.05 


2 


153 


06 


3.74 


156.19 


3.06 


160.05 


2.38 


164.01 


1.83 


168.78 


1 


39 


173.71 


1.09 


178.62 


0.87 


3 


193 


63 


4.61 


198.47 


3.69 


204.31 


2.77 


210.15 


2.08 


216.87 


1 


48 


223.37 


1.08 


229.48 


0.83 


4 


227 


51 


3.60 


231.93 


2.89 


237.25 


2.17 


242.58 


1.65 


248.81 


1 


22 


255.24 


0.98 


261.98 


0.84 


5 


262 


04 


4.58 


268.53 


3.67 


276.43 


2.78 


284.42 


2.11 


293.84 


1 


56 


303.36 


1.19 


312.62 


0.92 


6 


300 


74 


3.87 


306.90 


3.00 


314.08 


2.17 


321.07 


1.62 


329.02 


1 


15 


336.68 


0.85 


343.89 


0.66 


7 


332 


30 


3.48 


338.56 


2.81 


346.21 


2.16 


353.90 


1.62 


362.70 


1 


16 


371.30 


0.88 


379.95 


0.74 


8 


377 


94 


3.54 


384.60 


2.47 


391.56 


1.63 


398.06 


1.23 


406.14 


1 


02 


415.55 


0.92 


425.99 


0.80 


9 


401 


62 


2.99 


408.76 


2.81 


418.85 


2.48 


430.26 


2.06 


444.54 


1 


58 


458.90 


1.16 


471.83 


0.82 


10 


438 


27 


3.84 


447.38 


3.09 


458.37 


2.30 


469.24 


1.71 


481.02 


1 


12 


491.21 


0.75 


501.22 


0.69 


11 


471 


93 


3.67 


481.49 


3.05 


493.22 


2.31 


504.45 


1.60 


516.34 


1 


08 


528.34 


0.90 


541.75 


0.84 


12 


508 


90 


4.11 


519.69 


3.02 


530.98 


1.89 


540.66 


1.32 


552.67 


1 


14 


567.07 


1.01 


582.65 


0.88 


13 


548 


92 


2.53 


555.93 


1.88 


565.15 


1.75 


577.06 


1.69 


593.43 


1 


36 


609.46 


0.94 


623.80 


0.76 


14 


578 


30 


3.71 


590.30 


3.15 


605.61 


2.49 


621.19 


1.85 


637.61 


1 


14 


651.14 


0.77 


666.40 


0.86 


15 


607 


87 


3.18 


619.20 


2.90 


634.25 


2.36 


649.46 


1.72 


666.45 


1 


24 


684.86 


1.08 


705.77 


1.01 


16 


648 


21 


4.04 


662.06 


3.09 


676.56 


1.85 


687.85 


1.16 


701.46 


1 


05 


718.58 


0.96 


738.44 


0.96 


17 


689 


47 


3.71 


701.93 


2.49 


715.01 


1.79 


729.95 


1.67 


750.04 


\ 


ou 


768.60 


0.84 


787.10 


0.92 


18 


720 


75 


3.26 


733.19 


2.61 


749.99 


2.33 


769.03 


1.89 


790.64 


1 


24 


809.79 


0.90 


832.91 


1.06 


19 


752 


05 


2.99 


765.37 


2.82 


784.18 


2.43 


803.57 


1.75 


823.80 


1 


15 


845.00 


1.03 


871.16 


1.09 


20 


783 


56 


3.72 


800.55 


3.31 


821.32 


2.36 


839.23 


1.47 


859.00 


1 


21 


882.52 


1.05 


909.66 


1.10 


21 


825 


67 


3.79 


841.90 


2.79 


858.29 


1.68 


873.96 


1.48 


897.34 


1 


35 


921.88 


0.97 


948.35 


1.08 


22 


867 


99 


3.87 


883.99 


2.49 


900.69 


1.87 


920.13 


1.68 


943.76 


1 


14 


963.98 


0.79 


991.11 


1.14 


23 


896 


36 


3.53 


913.02 


2.83 


936.08 


2.56 


961.43 


1.96 


987.51 


1 


15 


1011.52 


0.96 


1042.90 


1.18 


24 


927 


94 


2.72 


943.81 


2.82 


967.15 


2.42 


990.17 


1.64 


1015.52 


1 


30 


1046.29 


1.18 


1081.47 


1.16 


25 


968 


77 


3.58 


988.78 


3.13 


1011.86 


2.00 


1030.38 


1.30 


1054.70 


1 


24 


1082.21 


0.94 


1114.22 


1.17 



NOTE - AU PjP values are in units of lO"^'^ s"^ 



The PjP values in our models exhibits a number of 
general trends. Firstly, we find that PjP generally de- 
creases with decreasing effective temperature, reflecting 
the diminishing effect of neutrino cooling and the increase 
in the cooling timescale due to the gradual decreasing of 
the thermal energy content. Second, trapped modes in 
the helium envelope have smaller PjP values since these 
modes are concentrated closer to the surface, where grav- 
itational contraction is still appreciable. As stated in §1, 
gravitational contraction acts shortening the pulsation pe- 
riods, causing trapped modes to have smaller P values. 
This trend is in agreement with previous studies on DB 
white dwarf pulsations (see, e.g., Bradley et al. 1993). 

Next, we shall briefly examine the effect of changing 
the various structural parameters on the rate of period 
change of our models. Factors that affect PjP in white 
dwarf models include the core composition, the surface 
chemical stratification, the stellar mass and the neutrino 
emission. Here we restrict ourselves to examining the ef- 
fects of the total stellar mass, the mass of the helium layer 
and the rate of neutrino emission. 

EfFect of the stellar mass: P/P in white dwarf stars is 
very sensitive mostly to the stellar mass. The effect of the 
stellar mass on the rate of period change has been explored 
by Kawaler et al. (1986) for DB white dwarfs, by Kawaler 
et al. (1985a) for DO white dwarfs, and by Bradley & 
Winget (1991), Bradley et al. (1993) and Bradley (1996) 
for DA white dwarfs. Lower mass models exhibit larger 



P/P values through the range of effective temperature 
of pulsating DB white dwarf stars (Kawaler et al. 1986). 
This can be understood on the basis that, if we consider 
a fixed Tes value, lower mass white dwarfs have larger 
luminosities and lower total heat capacity, and thus cool 
faster with a larger P. This can be also understood in 
terms of the simple, yet accurate enough for our purposes, 
cooling model of Mestel (1952). Within the framework of 
the Mestel (1952) cooling law, Kawaler et al. (1986) have 
derived a relation between P/P and the stellar mass: 

□ / \ 0.286 / 1,^ \ -1.190 

where /i is the mean molecular weight, /ig the molecular 
weight per electron, and A the atomic mass of the ions. In 
particular, it is convenient to note that from our numer- 
ical computations the magnitude of the period derivative 
varies by about a factor of 6 — 7 in the stellar mass inter- 
val 0.50 - 0.85Mq, as we show in Fig. H for DB models 
with a helium layer mass of « 8 x 10~^M». Each panel 
in the figure illustrates the situation at different effective 
temperatures covering the observed DBV instability strip. 
As we can see from the figure, the general theoretical ex- 
pectations from Eq. (0) arc borne out by our detailed nu- 
merical computations. In fact, the rate of period change 
for the less massive models is always larger than for the 
more massive ones. 
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Fig. 2. The P /P values versus the period, corresponding to DB white dwarf models with stellar masses from 0.50 to 
0.85 Mq (indicated with different line styles) and a helium content of Mhc = 8 x 10^'' M*. Left, centre and right panels 
display the results for effective temperatures of 30000, 25000, and 20000 K, respectively. Note: the strong minimum 
in P/P pointed with the arrow is not numerically inflicted (see text). 
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Periods and Relative Rates of Period Change — 1 modes) for selected models with M, = O.TOMq and Mhc = 9 x 10 ^M^, 







28075 K 


26997 K 


25992 K 


25055 K 


24066 K 


23033 K 


22077 K 


k 


P[s] 


P/P 


P[s] 


P/P 


P[s] 


P/P 


P[s] 


P/P 


P[s] 


P/P 


P[s] 


P/P 


P[s] 


P/P 


1 


125.06 


3.69 


128.96 


2.73 


132.89 


2.04 


136.88 


1.59 


141.50 


1.25 


146.84 


1.00 


152.29 


0.83 


2 


147.72 


2.93 


151.42 


2.25 


155.22 


1.69 


159.10 


1.33 


163.60 


1.04 


168.71 


0.82 


173.82 


0.67 


3 


189.01 


3.36 


194.33 


2.44 


199.56 


1.77 


204.62 


1.31 


210.11 


0.97 


216.09 


0.75 


222.06 


0.63 


4 


217.33 


2.67 


222.22 


1.98 


227.15 


1.49 


232.17 


1.19 


238.21 


0.99 


245.69 


0.86 


253.91 


0.77 


5 


254.45 


3.38 


261.71 


2.49 


269.05 


1.88 


276.45 


1.45 


284.90 


1.12 


294.10 


0.82 


302.39 


0.59 


6 


288.43 


2.60 


294.68 


1.89 


300.90 


1.41 


306.98 


1.05 


313.58 


0.78 


320.83 


0.62 


328.75 


0.60 


7 


318.29 


2.81 


325.80 


2.06 


333.05 


1.46 


339.93 


1.08 


347.65 


0.85 


356.89 


0.73 


366.92 


0.64 


8 


359.95 


1.90 


365.57 


1.38 


371.61 


1.18 


378.66 


1.08 


388.00 


0.95 


399.26 


0.77 


410.14 


0.59 


9 


384.26 


3.02 


394.63 


2.44 


405.74 


1.91 


417.08 


1.46 


429.36 


1.03 


441.26 


0.67 


451.90 


0.55 


10 


421.26 


2.73 


430.87 


1.98 


440.02 


1.37 


448.11 


0.91 


456.43 


0.70 


467.42 


0.72 


481.62 


0.74 


11 


452.58 


2.95 


463.50 


2.06 


473.38 


1.36 


482.64 


1.05 


494.02 


0.91 


507.86 


0.74 


521.77 


0.62 


12 


489.07 


2.35 


497.76 


1.49 


506.56 


1.28 


517.18 


1.19 


530.81 


0.99 


546.17 


0.75 


561.51 


0.67 


13 


518.80 


2.08 


529.51 


2.01 


542.42 


1.69 


555.59 


1.24 


568.74 


0.81 


581.56 


0.60 


596.47 


0.66 


14 


555.63 


2.93 


569.52 


2.19 


582.78 


1.48 


594.12 


0.96 


606.27 


0.80 


622.33 


0.75 


641.22 


0.75 


15 


581.77 


2.87 


595.88 


2.11 


609.22 


1.46 


622.58 


1.21 


639.63 


1.05 


659.93 


0.83 


681.23 


0.78 


16 


622.35 


2.44 


633.07 


1.35 


642.99 


1.14 


655.41 


1.12 


671.50 


0.91 


689.33 


0.70 


709.76 


0.77 


17 


657.81 


2.18 


671.06 


1.94 


686.81 


1.62 


702.30 


1.12 


716.70 


0.69 


731.67 


0.60 


752.84 


0.80 


18 


688.04 


2.70 


705.03 


2.24 


722.58 


1.62 


738.56 


1.12 


756.40 


0.93 


778.48 


0.78 


803.75 


0.84 


19 


719.30 


2.87 


737.04 


2.14 


753.25 


1.39 


768.50 


1.12 


788.46 


1.01 


812.59 


0.81 


840.19 


0.87 


20 


755.15 


2.95 


772.03 


1.77 


786.70 


1.31 


803.73 


1.23 


824.86 


0.95 


847.08 


0.71 


874.44 


0.88 


21 


789.76 


2.30 


804.40 


1.70 


821.99 


1.59 


841.50 


1.23 


861.54 


0.83 


883.89 


0.75 


913.96 


0.90 


22 


828.10 


2.13 


844.91 


1.95 


864.04 


1.49 


881.14 


0.98 


899.61 


0.83 


924.27 


0.76 


955.39 


0.90 


23 


857.68 


2.94 


880.47 


2.37 


902.55 


1.57 


921.73 


1.11 


944.67 


0.94 


970.05 


0.69 


1001.49 


0.92 


24 


888.14 


2.89 


909.34 


1.99 


927.99 


1.36 


948.99 


1.31 


976.68 


1.08 


1007.50 


0.83 


1044.05 


0.94 


25 


930.34 


2.70 


948.56 


1.56 


966.30 


1.39 


987.57 


1.17 


1009.91 


0.78 


1035.44 


0.76 


1072.25 


0.94 


26 


967.99 


2.51 


987.79 


1.90 


1010.83 


1.60 


1032.55 


1.04 


1054.27 


0.82 


1082.24 


0.73 


1119.22 


0.96 


27 


998.51 


2.41 


1021.85 


2.23 


1047.28 


1.59 


1069.80 


1.12 


1097.37 


0.99 


1128.35 


0.72 


1167.09 


0.97 



NOTE - All P/P values are in units of 10"^® s'^. 
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TABLE 7 

1 modes) for selected models with M, = 0.75Mq and Mho = 8 x lO^^M, 
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P 


si 

J 
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P [si 
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P fsl 
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P [si 


PIP 


P [si 
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1 


115 


07 


2.03 


118.26 


1.61 


121.98 


1.30 


125.83 


1.10 


130.99 


0.92 


135.77 


0.81 


141.33 


0.71 


2 


147 


18 


1.88 


150.91 


1.46 


155.16 


1.14 


159.35 


0.92 


164.61 


0.72 


169.11 


0.59 


173.95 


0.49 


3 


173 


02 


1.43 


176.28 


1.08 


179.89 


0.83 


183.44 


0.69 


188.17 


0.60 


192.78 


0.57 


198.63 


0.55 


4 


213 


57 


1.72 


218.49 


1.33 


224.12 


1.06 


229.91 


0.90 


237.66 


0.76 


244.68 


0.65 


252.41 


0.53 


5 


236 


16 


1.64 


241.58 


1.36 


248.16 


1.14 


255.04 


0.97 


264.13 


0.79 


272.26 


0.67 


281.49 


0.59 


6 


277 


38 


1.66 


283.33 


1.19 


289.30 


0.80 


294.32 


0.57 


300.19 


0.45 


305.85 


0.45 


313.36 


0.46 


7 


297 


53 


1.40 


303.14 


1.10 


309.88 


0.95 


317.46 


0.89 


328.60 


0.82 


339.35 


0.73 


351.84 


0.63 


8 


336 


46 


1.03 


341.60 


0.95 


348.47 


0.87 


355.90 


0.74 


365.24 


0.56 


372.78 


0.44 


380.97 


0.40 


9 


372 


43 


2.01 


382.53 


1.55 


393.38 


1.09 


402.32 


0.70 


411.34 


0.48 


419.17 


0.44 


429.07 


0.45 


10 


396 


26 


1.41 


403.35 


0.99 


410.92 


0.79 


419.88 


0.85 


435.03 


0.87 


450.46 


0.79 


468.69 


0.70 


11 


435 


31 


1.27 


442.23 


0.91 


450.48 


0.81 


459.61 


0.71 


471.35 


0.55 


481.21 


0.46 


493.41 


0.51 


12 


457 


41 


1.47 


467.69 


1.40 


480.84 


1.17 


494.06 


0.93 


509.82 


0.68 


522.69 


0.55 


538.12 


0.58 


13 


495 


52 


1.50 


505.36 


1.11 


515.23 


0.75 


524.02 


0.60 


536.72 


0.60 


550.82 


0.62 


570.27 


0.68 


14 


528 


02 


1.73 


539.38 


1.17 


550.72 


0.86 


562.24 


0.74 


577.67 


0.61 


591.10 


0.51 


608.65 


0.63 


15 


556 


80 


1.27 


566.46 


1.06 


579.67 


1.03 


594.65 


0.91 


613.88 


0.69 


629.75 


0.56 


650.51 


0.69 


16 


583 


05 


1.41 


595.15 


1.26 


609.90 


1.02 


624.25 


0.79 


642.26 


0.66 


660.27 


0.66 


685.76 


0.76 


17 


619 


08 


1.47 


631.67 


1.16 


644.78 


0.80 


656.76 


0.66 


673.88 


0.62 


690.99 


0.58 


715.23 


0.74 


18 


659 


60 


1.62 


672.23 


1.00 


684.72 


0.78 


698.55 


0.73 


717.55 


0.60 


733.91 


0.51 


759.14 


0.79 


19 


683 


63 


1.48 


697.58 


1.24 


715.78 


1.12 


734.86 


0.90 


757.29 


0.64 


776.21 


0.57 


804.98 


0.81 


20 


715 


69 


1.25 


729.03 


1.16 


746.06 


0.96 


762.71 


0.77 


785.46 


0.72 


808.91 


0.69 


841.61 


0.82 


21 


746 


97 


1.59 


763.68 


1.28 


781.17 


0.88 


797.14 


0.72 


819.76 


0.66 


841.48 


0.60 


874.08 


0.84 


22 


782 


59 


1.56 


797.99 


1.07 


814.10 


0.86 


832.45 


0.82 


857.35 


0.65 


878.77 


0.56 


912.26 


0.85 


23 


818 


75 


1.44 


833.23 


1.01 


851.15 


0.93 


870.18 


0.76 


893.50 


0.60 


915.85 


0.60 


951.68 


0.85 


24 


846 


51 


1.55 


865.32 


1.35 


888.22 


1.06 


908.75 


0.75 


933.41 


0.62 


956.87 


0.57 


994.35 


0.89 


25 


873 


67 


1.35 


891.63 


1.25 


913.03 


0.97 


935.22 


0.90 


968.32 


0.80 


997.63 


0.65 


1038.55 


0.90 


26 


915 


25 


1.52 


933.00 


1.04 


950.10 


0.76 


969.13 


0.74 


995.97 


0.63 


1022.32 


0.65 


1065.09 


0.86 


27 


952 


28 


1.64 


971.17 


1.09 


992.43 


0.95 


1014.28 


0.73 


1039.24 


0.55 


1063.09 


0.55 


1105.62 


0.92 


28 


980 


15 


1.47 


1000.22 


1.27 


1026.48 


1.08 


1051.18 


0.80 


1082.11 


0.67 


1109.78 


0.57 


1154.29 


0.93 
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Fig. 3. The P/P values versus the period, corresponding to DB white dwarf models with a stellar mass of 0.60 Mq 
and a content of helium of 1 x 10"'' (solid lines), 9 x 10""' M* (dashed lines), and 8 x 10"^M, (dotted lines). Left, 
centre and right panels display the results for effective temperatures of 30000, 25000, and 20000 K, respectively. P/P 
values corresponding to models with thick helium envelopes are only marginally smaller as compared with the case of 
thin helium envelopes. 
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Fig. 4. The evolution of the relative rate of period change 
including (continuous lines) and ignoring (short-dashed 
lines) the effect of neutrino, in terms of the effective tem- 
perature, corresponding to a sequence of O.GO-Mq white 
dwarf models (Mho = 8 x IQ-^M*). 



Fig. 121 further emphasizes the complex structure char- 
acterizing the P/P-distribution: no matter the stellar 
mass is, the diagrams show abundant local maxima and 
minima. The mode with fc = 14 corresponding to the 0.50- 
Mq white dwarf model at T^s ~ 30000 K (pointed with 
the arrow in left panel of the figure) is striking. The very 
pronounced minima in the rate of period change of this 
mode is not the result of a numerical glitch; we have ver- 
ified the reality of this by examining their oscillation ki- 
netic energy and eigenfunctions, and we have found that 
this mode is strongly trapped in the helium-rich envelope. 
As we have already anticipated, trapped modes like this 
one must have a small P/P value, reflecting primarily the 
contraction of the outer layers, rather than carrying in- 
formation about the processes linked to the cooling of the 
white dwarf (Bradley et al. 1993). We have verified that 
almost all the trapped modes in the helium-rich envelope 
of our models correspond to local minima in P/P. 
Effect of tlie fielium layer mass: Now we examine the 
possible dependence of the P/P values on the thick- 
ness of the helium layer, holding the stellar mass fixed. 
We have considered the following values for the helium 
mass: 1 x 10"'' M,, 9 x 10^'' A/,, and 8 x IQ-^M^. It is 
worth recalling that in the sequence of DB models with 
Mho ~ 1 x 10^'^M, the envelope is characterized by an 
initially double-layered configuration which evolves to a 
single-layered structure before the model reaches the red 
edge of the DBV instability strip (Althaus & Corsico 
2004). The effects of varying the mass of the helium layer. 



keeping the total mass constant, are illustrated in Fig. 13 
This figure shows that the P/P values corresponding to 
models with massive helium envelopes are slightly smaller 
than the case of thinner ones. This effect can be explained 
on the basis that our models with massive helium layers 
cool slightly more slowly as compared with models charac- 
terized by thinner helium envelopes. This is understood on 
the basis that models with massive (and also more trans- 
parent) helium envelope are characterized by lower central 
temperatures which in turn, implies that such models have 
initially an excess of internal energy to get rid of — see 
Tassoul et al. (1990). This leads to longer cooling age at 
this stage and thus smaller rate of period variations. 
Effect of the neutrino emission: Finally, we shall examine 
the effect of the neutrino emission on the P/P values. We 
begin by examining Fig.^ in which we show the evolution 
of the rate of period change in terms of the effective tem- 
perature corresponding to a sequence of O.6O-Af0 white 
dwarf models with a helium content of Mne = 8 x 10~^M». 
Continuous lines correspond to the case in which neutrino 
emission has been included in our DB models, and short- 
dashed lines correspond to the situation in which neutrino 
emission has been ignored^. At the hot edge of the DB 
instability strip, the values of P/P considering neutrino 
emission are greater by a factor of about 5 as compared 
with calculations in which neutrino emission has not been 
taken into account. At Teg « 25000 K the factor is of 
about 2.5. In Fig. we illustrate P/P versus periods at 
three effective temperatures, for the case in which neu- 
trino losses have been considered (continuous lines) and 
for the situation in which neutrino emission has been ig- 
nored (dashed lines). Clearly, the P/P values are very 
sensitive to neutrino emission, particularly at the high ef- 
fective temperatures characterizing the blue edge of the 
DBV instability strip. Thus, we essentially recover the re- 
sults reported by Winget et al. (2004). We conclude, in 
agreement with Winget et al. (2004) 's claims, that at high 
effective temperatures within the DB instability strip the 
eventual detection of the rate of period change in DBV 
white dwarfs could allow the astronomers to constrain the 
production rates of plasmon neutrino emission. 

4. Summary 

The rate of period change is the most exciting observable 
quantity in pulsating white dwarf stars because it poten- 
tially provides a direct measure of the rate of cooling of 
white dwarfs, giving astronomers the possibility of infer- 
ring the age of the galactic disk in the solar neighborhood 
(Winget et al. 1987). Also, measurement of P in white 
dwarfs offers an unique opportunity to place interesting 
constraints on neutrino physics (O'Brien & Kawaler 2000; 
Winget et al. 2004) and to test any additional sink of 

^ It is worth mentioning that we have turned off neutrino 
emission at a high effective temperature (~ 60000 K). As a re- 
sult, the model has adequately readjusted its thermal structure 
long before the DBV regime is reached. 
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1 


81 


130.85 


1.40 
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1.15 
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0.94 


145.06 


0.79 
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0.65 


2 
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97 
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146.53 


1 


53 


150.59 


1.18 


154.64 


0.95 


159.26 


0.76 


163.89 


0.64 


169.15 


0.55 


3 


183 


88 


2.04 


188.41 


1 


48 


193.25 


1.07 


197.86 


0.84 


203.16 


0.70 


208.78 


0.63 


215.67 


0.58 


4 
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81 


1.77 


212.37 


1 


37 


217.80 


1.12 


223.67 


0.99 


231.05 


0.87 
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0.77 


248.49 


0.66 


5 


246 
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253.53 
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1.26 


268.46 


0.97 


276.19 


0.69 


282.82 


0.49 


289.35 


0.38 


6 
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1.67 


281.04 


1 


23 


287.05 


0.89 


292.70 


0.70 


299.52 


0.64 


307.70 


0.65 


318.76 


0.64 


7 
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1.78 


312.76 


1 


24 


319.61 


0.94 


326.72 


0.82 


335.56 


0.72 


344.83 


0.61 


355.10 


0.50 
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44 


1.33 


346.67 
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355.13 


1.10 


364.26 


0.91 


374.49 


0.70 


383.97 


0.54 


394.41 


0.48 
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00 


2.22 


381.23 
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69 


392.44 


1.20 


402.31 


0.83 


412.13 


0.61 


422.03 


0.55 


434.67 


0.55 


10 
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0.74 


425.52 


0.74 


437.16 


0.77 
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Fig. 5. The relative rate of period change in terms of the period corresponding to the same models as in Fig. 0] at 
effective temperatures of « 30000, « 25000 and ~ 20000 K. Continuous (dashed) lines correspond to computations 
including (ignoring) the effect of neutrinos. 



energy (Isern et al. 1992; Corsico et al. 2001a). In addi- 
tion, the small magnitude of P characterizing pulsating 
white dwarfs implicates that its detection would also im- 
pose strong constraints on the presence of other mecha- 
nisms acting on shorter time scales and that can modify 



the pulsation periods, such as stellar rotation (Kawaler et 
al. 1985b) and binary orbital motion (Kepler et al. 1991). 

Rates of period change have been detected for two 
members of the ZZ Ceti variable class (G117-B15A and 
R548) and only for a single DOV white dwarf (PC 1159- 



12 Corsico & Althaus: The rate of period change in DBV white dwarf stars 

TABLE 9 



Periods and Relative Rates of Period Change {£ — 1 modes) for selected models with M» — 0.85Mq and A/hc = 8 x 10 ^Af, 
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NOTE - All P/P values are in units of 10"^'' 8"^ 



035). In the case of ZZ Ceti stars, the observed values for 
the rate of period change are in good agreement with the 
theoretical expectations. As for PG 1159-035, however, 
the observed value is largely in excess as compared to the 
theoretical predictions (see Costa et al. 1999). For DBV 
stars no detection has yet been made, but this situation 
could soon be changed by the ongoing observations of the 
hot DBV EC 20058 by Sullivan et al. (2004). In a recent 
paper, Wingct ct al. (2004) have called attention to the 
potential of employing DBV white dwarf stars to place 
constraints on the plasmon neutrino emissivity. Prompted 
by Wingct ct al. (2004) 's claims, wc have computed in 
this work the change of period rate in pulsating DB white 
dwarf stars employing the new stellar models recently pre- 
sented in Althaus & Corsico (2004). These stellar models 
incorporate an updated input physics in addition to a self- 
consistent treatment of the time-dependent element diffu- 
sion to assess the shape of the chemical interfaces — an 
issue of fundamental importance in the context of white 
dwarf pulsations. In view of the lack of modern tabula- 
tions of P/P in DBV stars in the literature, we have in- 
cluded tables providing P and P/P values corresponding 
to dipole modes (^ = 1) for different stellar masses, effec- 



tive temperatures and thickness of the helium envelope. 
We have also examined the effects of varying the neutrino 
emission rate on the P/P values. Wc find, in agreement 
with previous studies, that the rate of period change is a 
function mainly of the total stellar mass — there is more 
than a factor 6 between the P/P values for 0.50- and 
O.85-M0 models. At variance, the rates of period change 
are rather insensitive to the thickness of the helium enve- 
lope. We also find that, within the DBV instability strip, 
P/P is strongly dependent on the neutrino energy losses. 
Specifically, the P/P values are considerably smaller (by 
a factor of about 3 — 5) when neutrino losses are ignored in 
our models. Our results are consistent with the conclusion 
arrived at recently by Winget et al. (2004) about the fea- 
sibility of quantitatively constraining the production rate 
of plasmon neutrinos by employing pulsating DB white 
dwarfs. 
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